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Abstract 

A chiral fourth generation is a simple and well motivated extension of the standard 

model, and has important consequences for Higgs phenomenology. Here we consider a 
scenario where the fourth generation neutrinos are long lived and have both a Dirac 
and Majorana mass term. Such neutrinos can be as light as 40 GeV and can be 
the dominant decay mode of the Higgs boson for Higgs masses below the W-boson 
threshold. We study the effect of the Majorana mass term on the Higgs branching 
fractions and reevaluate the Tevatron constraints on the Higgs mass. We discuss the 
prospects for the LHC to detect the semi-invisible Higgs decays into fourth generation 
neutrino pairs. Under the assumption that the lightest fourth generation neutrino is 
stable, its thermal relic density can be up to 20% of the observed dark matter density 
in the universe. This is in agreement with current constraints on the spin dependent 
neutrino-neutron cross section, but can be probed by the next generation of dark matter 
direct detection experiments. 



1 Introduction 



The fermionic matter content of the standard model of particle physics is organized into 
three families with chiral couplings to the electroweak SU(2)xU(l) gauge bosons. The 
minimum of three generations is required by the CKM mechanism for CP violation in the 
standard model, and has of course been confirmed in the last decades by the discovery of 
the top quark, the last missing member of the third family of quarks and leptons. On the 
other hand, there is no immediate upper limit on the number of fermion generations, such 
that a fourth generation is a well motivated extension of the standard model that should 
be searched for by experiments. 

Phenomenology of a fourth generation of quarks has been studied exhaustivel}j^ Direct 
searches for additional up and down-type quarks by the Tevatron and very recently at the 
LHC constrain their masses to be above 335 GeV and 372 GeV respectively |3lllj. 

These searches however depend on the decay modes of the fourth generation quarks [5] . 
A less model dependent way to constrain fourth generation quarks is through their effects 
on electroweak precision constraints [HI d EJ [2] and on the Higgs production cross section 
O [9l [To] . The latter has been used recently by the Tevatron to exclude Higgs masses from 
131 GeV-204 GeV in the presence of a fourth generation |llj . 

Precise measurements of the Z-boson width at LEP excludes the existence of a fourth 
neutrino flavor with a mass below M^/2. It follows that the fourth generation neutrino 
must have a Dirac or a Majorana mass term. In addition the LEP experiments have 
searched for fourth generation leptons that decay to SM leptons, and also for charged stable 
leptons, roughly constraining them to be heavier than 100 GeV. For a recent analysis of 
these constraints for mixed Dirac Majorana neutrino masses see [l^. Fourth generation 
Majorana neutrinos are also constrained by limits on lepton flavor violation, assuming a 
nonzero mixing with the first three generations of fermions |13j . 

If the mixing of fourth generation fermions with the standard model is very small, 
respectively the lightest quark and lepton become long-lived |14j. It was recently pointed 
out that such long-lived quarks not only have distinct phenomenological signatures, but 
also can assist electroweak baryogenesis scenarios by protecting baryon number from being 
washed out by sphaleron processes [T^]. In this case also the LEP constraints on fourth 
generation neutrinos must be revisited |16j . 

A long lived fourth generation neutrino can have important consequences for Higgs 
phenomenology, since the coupling to the Higgs boson is proportional to the neutrino 
mass. For niu < mfi/2 the decay into neutrino pairs easily outweighs the Higgs to bb width 
and can be the dominant Higgs decay mode for Higgs masses below 2M]y- In section 2 we 
review the properties of mixed Dirac-Majorana neutrinos, using a two component spinor 
formalism, and derive the couplings to the Higgs and Z-bosons. In section 3 the effects on 
Higgs production and decay are discussed, paying attention particularly to the effects of 

^See e.g. jl] for a review, an extensive list of references can also be found in [2]. 
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the Majorana mass term that splits the two neutrino states. We proceed by discussing the 
colhder phenomenology of Higgs bosons in this scenario in section 4. Finally we calculate 
the relic density and direct detection constraints on a stable fourth generation Majorana 
neutrino in section 5, before presenting our conclusions in section 6. In the appendix we 
translate the results of section 2 to four component spinors. 



2 Fourth Generation Neutrino Masses and Couplings 

It is instructive to first review the masses and couplings of mixed Dirac-Majorana neutrinos, 
see also |17j for a similar analysis. We work in the two component spinor formalism, 
following the conventions of [18]. The Higgs doublet and the fourth generation neutrinos 
and charged leptons are denoted by 

^=(^^), Vu, Vi, (1) 

with hypercharges {Y = Q—T3) assigned as +5, — 0, +1 respectiveljj^ . Gauge invariance 
further implies that the charged Weyl fields X£, Ve can be combined into a Dirac spinor 

The two types of SU(2) singlets that can be constructed are (j)^L and (f^rL, with 
hypercharges —1 and 0, and r = —ia2- Therefore the allowed Yukawa couplings are 

£ D (t^^L 7]e + Yn (frL ri^ + h.c. (2) 

These are the usual charged and neutral lepton Yukawa couplings. Note that we have 
neglected any mixing with the first three generations of leptons. The Higgs field, in unitary 
gauge, is decomposed as 

1 



, ^ + , (3) 

with V = y/2mw/g » 174 GeV. 

Here we are mainly interested in the interactions of the physical Higgs boson with the 
neutrinos. Dropping the subscripts on Xv and r/i,, the relevant part of the Lagrangian 
becomes 

+ -^l^ JXV+ yn4D + -^'T' J X'V + [VV + V V j , (4) 

where we have introduced the Majorana mass term, and defined m^D = YnV. The resulting 
mass matrix can be diagonalized using the Takagi decomposition [18] . Note that in general 



^Note that all undaggered fields are left-handed, i.e. rii is the left-handed positron. 
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the masses are not the eigenvalues of the mass matrix Ai , but the positive square roots of 
the eigenvalues of Ai'^M. The masses and corresponding eigenstates are given by 



Mi,2 = \l^+ml^±-, (5) 

^^1 = i {ceX - sev) , (6) 
^^2 = sex + cer] , (7) 



with 



Ma m4D V Ma ^ ^ 

The physically relevant quantities are the two Majorana masses Mi^2; which we take as the 
two independent free parameters of the model. The coupling to the Higgs boson, which is 
proportional to m^D, can be rewritten in terms of Mi^2 as follows: 

m4D = IViMi + M2Y - (M2 - Mi)2 = VM1M2 = Mi^l^^. (9) 

Note that for fixed Mi, which is bounded by LEP to be larger than about 40 GeV, it 
increases with the square root of the ratio of the masses. To remain in the perturbative 
regime, we require m^D ^ 400 GeV, which implies that M2 < 4 TeV. 

The couplings to the Higgs boson are now obtained using equations ^ and ([T]): 

C = ixv + h.c.) = J^/i {2ceseNf + 2ceseNl + 2i{sj - cDn^N^ + h.c.) . (10) 

We have extracted factor of 1/2 to obtain the canonical normalization for the Majorana- 
Higgs couplings. In terms of Mi ^2; the coupling of the lightest neutrino to the Higgs is 
given by 

m4D„ Ml 2M2 

In the limit M2 — )■ Mi, the sum of the partial widths should agree with the width of the 
Higgs decaying into a single Dirac neutrino. The sum of the widths is proportional to 

\ [Aclsl + Aclsl) + (4 - slf = (4 + .2)2 = 1 , (12) 

so it is independent of the mixing angle, in the limit where the phase space factors can be 
neglected, i.e. for mh ^ 2Af2. The factor 1/2 in front of the first term is from the phase 
space for identical particles. The 1/2 in the coupling (10) is compensated for the NiNi 



and N2N2 channels since each field can be contracted with each final state. For the N1N2 
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channel the factor 1/2 cancels the factor of two in the coupling, such that the relative 
factors emerge as shown above. 

For completeness, we also give the couplings to the Z boson. In two component notation, 
they are obtained from 



= ^ 



(13) 



In four component notation, this translates into axial couplings of the Z boson to NiNi 
and N2N2 and a vector coupling to N1N2, as discussed in more detail in the appendix. 

Before moving to the next section, let us briefly comment on the phenomenology of 
the fourth generation charged leptons in this model. Direct searches at LEP constrain the 
mass of the charged lepton to be larger than 100.8 GeV [23j, while electroweak precision 
tests constrain the mass difference between the charged and neutral lepton, \m£_^ — m,/^\ < 
140 GeV [2J, for a Dirac 1/4. 

Pair production of £4^4 pairs with subsequent decays ^4 — ?■ WNi^2 lead to signatures 
with leptons, jets and missing energy. In the presence of a nonzero Majorana mass term, 
same sign lepton pairs can be produced in the process 

pp ^ if if i^i^ + ^ + X, (14) 

where X denotes the proton remnants from the W fusion process. In addition charged- 
neutral pairs i4:Ni^2 can be produced and lead to signals with one charged lepton and 
missing energy. Some of these processes were considered e.g. in [T21[IZ]. For a recent study 
of same sign leptons see also |19) . however in our case the stable A^i and the neutrinos from 
W decays lead to much larger missing energy. 



3 Higgs Production and Decay 

Higgs production through gluon fusion is enhanced in the presence of a heavy fourth gen- 
eration. For TJifi < 200 GeV, the enhancement is well approximated by multiplying the 
standard model production rate by a factor = 9 [3191 [10], and is largely independent of 
the structure of the lepton sector. If Higgs decays to the fourth generation are not allowed 
kinematically, this greatly strengthens the bounds on the Higgs mass from the Tevatron 
experiments, excluding the 131 — 204 GeV region |llj . 

The light Higgs mass range can be resurrected in the presence of a long lived fourth 
generation due to the invisible decays h — )• 1^41^4. The tree level decay width into a fourth 
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Figure 1: Higgs branching fractions as function of the Higgs mass. Sohd hnes correspond 
to a fourth generation Dirac neutrino with mass Mi = 45 GeV, dashed lines are the 
corresponding standard model branching fractions. Above the h — )■ WW threshold the 
branching fractions converge to their corresponding standard model values, except for the 
gluon channel that remains enhanced. 



generation Dirac neutrino with mass Mi is given by 



3/2 



r,, = T{h ^ v^v^) = Irrih 1 f . 15 

For a light Higgs this competes with the decay into b quarks, and the ratio of the width is 
approximately given by 

The decay into fourth generation neutrinos therefore dominates the total width of the 
Higgs. In Figure Q we show the branching fractions of the Higgs boson in the presence of 
a purely Dirac neutrino z^4 with a mass of 45 GeV. We have also included the enhancement 
of the h ^ gg channel due to the fourth generation quarks. Our findings for the case of 
a Dirac neutrino agree with the results of previous studies |20^ I21j . The SM branching 
fractions of the Higgs were obtained using HDECAY [22] . 

The situation becomes more complicated in the presence of a nonzero Major ana mass. 
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In this case, we find that the partial widths are given by 




r(. N,N,) = ^'^^Jfifl (i - ^^^^1 ,/a (* . ^1 . (19) 

167r (M2 + Mi)^ \ J \ m-h. 

The definition for the triangle function A can be found in the appendix. As discussed above, 
in the Dirac limit, where M2 — )• Mi, the width into N1N2 goes to zero, and Tn +T22 = ^uu- 
On the other hand, in the limit where M2 ^» Mi,mh, only the decay into NiNi is possible. 
It is interesting to note that in this case the factor 1/2 in the width is overcompensated by 
the factor {2M2f/{Mi + M2f 4, and we have that 

Tih ^ iViiVi)|,,^>>^,^ = 2T{h ^ (20) 

This can also be understood in terms of an effective theory. Integrating out the right- 
handed neutrino, the Majorana mass for the left-handed neutrino comes from a term 

^^(0L)(,^L) = i^(0L)(0L) ^ Im,N,N, + ^-^^hN,m , (21) 

i.e. the coupling to the Higgs is enhanced by a factor of two compared to the case of a 
Dirac mass term. 

We can therefore conclude that the presence of a Majorana mass term in the fourth 
generation neutrino sector can increase the invisible width of the Higgs boson by up to a 
factor of two. It is decreased compared to the Dirac case in a narrow region of parameter 
space, where M2 is just large enough to forbid Higgs decays into N2N2. The suppression 
is 50% at most, i.e. Fn > T,yy/2. 

To illustrate the effects of a nonzero Majorana mass, in figure [2] we show the branching 
fractions of the Higgs boson for two characteristic cases. We always take Mi = 45 GeV to 
simplify the comparison with the Dirac case. For the left plot we take M2 = Mi + 10 GeV 
such that all decay channels are open. The right plot is valid in the limit where 
decouples, i.e. M2 S> Afi. 

Comparing with the Dirac case, one can see that the suppression of the standard model 
Higgs decays is somewhat weaker in the presence of a small Majorana mass, while it is 
stronger for the case of a large Majorana mass term. The N2N2 channel suffers from 
phase space suppression compared to the NiNi channel, but will be comparable whenever 
nih » 2M2. The A^iA''2 channel is always suppressed, either by the (M2 — Mi)^ factor in the 
width, or by phase space when M2 3> Mi. For moderate mass splittings M2 — Mi ~ 20 GeV 
the branching fraction can reach up to 8%, for Higgs masses below the WW threshold. 



7 



100 



1.000 
0.500 

0.100 
B 0.050 

^ 0.010 
0.005 

0.001 



150 200 
m,,[GeV] 



300 



100 



150 200 
m;,[GeV] 



300 



Figure 2: Higgs branching fractions in the presence of stable fourth generation Majorana 
neutrinos. In the left plot we assume a small mass splitting, Mi = 45 GeV and M2 = 
55 GeV. The right plot shows the limit where decouples from the theory. 
Thick red lines are the Higgs branching fractions into NiNi (solid), N2N2 (dashed) and 
A''iA''2 (dotted). The black lines are the WW (solid) and ZZ (dashed) channels, while the 
blue (light grey) lines are the GG (solid) and bb (dashed) channels. Decays into charm, 
strange and tau pairs are not shown, other decays have branchings smaller than 5 x 10~^. 



4 Higgs Phenomenology with mixed A^i, A^2 

The limit on the neutrino mass of Mi > 45 GeV only applies for stable Dirac neutrinos. 
For a pure Majorana A^i instead one obtains Mi > 39.6 GeV [23]. In the mixed case 
the constraints on Mi can be further reduced, since the coupling to the Z boson can be 
suppressed by the mixing angles. In addition one has to take into account LEP constraints 
on N2 ^ Ni + decays. This was recently analyzed in [16]. The allowed regions for Mi, 
M2 are roughly given by: 

I Ml > 30 GeV M2 - Ml < 20 GeV , (22) 

II Ml > 40 GeV M2 + Ml > 170 GeV . (23) 

For previous studies of Majorana neutrino production through Higgs boson exchange, see 
also i24i|25j and tlTj. 

4.1 Higgs Mass Bounds 

In the presence of a fourth generation, the production of Higgs bosons in gluon fusion is 
enhanced by roughly a factor of 9 [7J. On the other hand the associated production of a 
Higgs with a Z or W boson is not significantly changed. 
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Figure 3: Suppression of the branching fraction Br(h — ?■ WW*) relative to the SM branching 
fraction, as a function of rrih and Mi, for M2 = 4000 GeV. Shown are the contours for 
the ratio Br4G(h — )■ WW*)/BrsM(h — )• WW*). Below the red contour, the suppression 
is stronger than the enhancement of Higgs production in gluon fusion. The suppression 
in the region where the decay h — )• NiNi is not allowed is due to the fourth generation 
enhancement of the Higgs decay to gluon pairs. 



The LEP experiments pose stringent bounds on the invisibly decaying Higgs boson, 
such that nih > 114 GeV continues to hold in any scenario considered here. In the low 
mass region, mh < 130 GeV, the strongest constraints come from ZH and WH channels. 
These limits are not enhanced by the increased gluon fusion rate, such that the low mass 
region remains viable. 

Above m/i = 130 GeV, the pp ^ h ^ WW* mode is strongly enhanced and, in the 
absence of new Higgs decay channels, excludes a Higgs in the mass range 131 GeV < ruh < 
204 GeV [11] . The suppression of the WW width from Higgs decays into fourth generation 
neutrinos is crucial to resurrect Higgs masses between 131 GeV and 2M\y. The relevant 
condition is 

^ ^ Bvsuih^WW*) - ^^■^^^'"'^^^ ^^^^ 

where (Jh,exduded is the Higgs production cross section excluded by the combined Tevatron 
limit, normalized to the standard model production cross section. The suppression of the 
h — )■ WW branching fraction for different values of Mi is shown in Fig. |3] for the case when 
M2 decouples. For Mi = 50 GeV, we obtain an upper bound on the light Higgs mass of 

rrih < 156 GeV , (25) 
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which is vahd for M2 ^ Mi, i.e. for a large Majorana mass term. This result is quite 
stable under small variations of Mi . Lowering the mass decreases the coupling to the Higgs 
but enlarges the available phase space, while increasing the mass has the opposite effect. 
This can also be seen in Fig. [Sj where the contours of equal suppression are approximately 
straight vertical lines, away from the threshold for h ^ NiNi. In the regime where 
M2 — Ml <^ M2, i.e. for a small Majorana mass, the h — )• WW* width is slightly less 
suppressed, such that Higgs masses down to 152 GeV become excluded. 

In the large Higgs mass region, the branching into on-shell WW pairs is only slightly 
reduced by the neutrino decay modes, such that the lower bound for the high mass region 
is not significantly changed from > 204 GeV obtained in |llj . 

In region I the decays h — t- N1N2 and h — t- N2N2 are possible, and offer a possibility 
to discover the otherwise invisible decays of the Higgs bosons through the decays of N2 to 
A'^i, through an off-shell Z or Higgs boson. 



4.2 N2 Decays 

In addition to the decays — ^ NiZ^*'^ that have been studied in the literature, one 
must also take into account decays of into on- or off-shell Higgs bosons. The relevant 
couplings to the Z and Higgs are 

g Mz 2M2 

= vrvWTM-2 

-^i4-sl)-V-2^^ (27) 



^/2v " " V M2 + M1 

When the mass splitting is small, the coupling to the Higgs is suppressed. In addition, 
since both decays are into virtual particles here, the small width of a light Higgs boson 
makes this channel even more suppressed. When M2 is large enough such that both decays 
can happen on-shell, we will roughly have that Th/Tz ~ MI/M"^. In this regime the Higgs 
will mostly decay invisibly. We can therefore conclude that, as far as Higgs physics is 
concerned, it suffices to consider the decays NiZ. 

For direct A''2 searches at the LHC, other parameter regions might become important. 
In particular for M2 — Mi > rrih and either rrih ~ 155 or ruh > 200, we will have a situation 
where the Higgs decay to W~^W~ is significant. This can lead to the following interesting 
signals: 

pp ^NiN2 ^W+W- (28) 
pp ^N2N2 ^ W+W- Z + (29) 

These channels might offer improved sensitivity for LHC searches in the high mass regime. 
With production cross sections in the (10 — 100) fb range at the ^/s = 7 TeV LHC it will 
still be a very challenging search. 
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4.3 Higgs Phenomenology with Small M2 



We first consider the light Higgs case. Since M2 — Mi < 20 GeV, both h — )• N1N2 and 
h — )• N2N2 are kinematicahy accessible. N2 decays to A^i via an off-shell Z boson, while 
decays mediated by an off-shell Higgs are suppressed. The most promising channel therefore 
is 



One can look for the Higgs boson in the di-lepton and four lepton channel with missing 
energy. For Higgs bosons produced in the gluon fusion mode, the leptons will have very 
small transverse momenta, and fail most experimental cuts. It is more promising to consider 
associated Zh and Wh production. Not only will the Z and W bosons provide a trigger, 
we can also expect a moderate boost for the Higgs, such that the final state leptons pass 
some mild pT cuts. 

The production cross section for Zh and Wh at the Tevatron, for a light Higgs, are 
0(100 fb). Channels with the Higgs going to four leptons are suppressed too much by the 
Z boson branching fractions. The most promising channel is 



with a cross section of 0(1 fb) at the Tevatron. The most recent measurements from the 
DO and CDF experiments constrain the cross section for pp ^ 3£ + to be smaller than 
0.1 pb, still two orders of magnitude above the expected signal. 

At the LHC with ^/s = 7 TeV the WH and ZH production cross sections can reach 
up to 0.5 pb. Backgrounds for this channel come from WZ and Wj* production as well 
as from tt production, where additional leptons are produced in decays of heavy flavors. 
A recent CMS study [27j shows that, with a moderate cut on missing transverse energy, 

> 50 GeV, and requiring that the invariant mass of same flavor lepton pairs is away 
from the Z pole, the backgrounds can be reduced to the 30 fb level, comparable to the 
expected signal rates. 

Finally the conventional Higgs search channels remain viable, although modified by a 
factor 9 x r4G'/rsM. In the case that the Higgs is found in one of those channels, the 
measured production cross section will be different from the standard model prediction. It 
is then crucial to look for the semi-invisible Higgs decays to understand the origin of this 
deviation from standard model predictions. 

4.4 Higgs Phenomenology with Large M2 

Here we assume that we are in case (II), but without decoupling A''2. In the small rrih region, 
the Higgs decays only to NiNi, which is invisible. Finding the Higgs in this regime will be 
difficult, since the conventional search channels are significantly suppressed, in particular 



N2 i+r + $ . 



(30) 




(31) 
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those based on associated Higgs production. The invisibly decaying Higgs may however be 
observed in the weak boson fusion channel at the LHC |28j . 

For the region 200 GeV < ruh < 300 GeV, the Higgs decays into pairs of electroweak 
gauge bosons dominate. The N1N2 and N2N2 channels suffer from phase space suppression 
in this regime, such that their branching fractions are only at the 1% level. 

The now decays into an on-shell Z boson. The relevant leptonic channels are 



The di-lepton channel receives contributions both from N1N2 and N2N2 decays, when one 
of the Z bosons decays invisibly. The dominant SM backgrounds for this channel are 
ZZ and WW pair production. In addition, there is also an irreducible background from 
electroweak N1N2 and N2N2 production, mediated by a Z boson. The signal rate for this 
channel at the 7 TeV LHC is about 



while it is largely invisible at the Tevatron. At the LHC the background from W^W~ pro- 
duction can be reduced by requiring the invariant mass of the lepton pair to be close to mz ■ 
However the background from ZZ — t- i~^£~ + ^ is irreducible and, with about 170 fb [|29j, 
a few times larger than the potential signal. This background can be measured accurately 
at the LHC using ZZ — t- 4^ decays. Direct N1N2 production contributes 0(20 fb) to the 
£+£- + ^ background. With 5 fb'^ at the 7 TeV LHC a significant excess in this channel 
might be observable, depending on the precise values of and Mi^2- 

5 Fourth Generation Neutrino Dark Matter 

For the above analysis to be valid, it is sufficient if the A'^i lifetime is long enough to escape 
the LEP detectors before they decay. However it is also conceivable that A^i is stable on 
cosmological timescales, and thus contributes to the dark matter of the universe. 

It is well known that stable fourth generation Dirac neutrinos make very bad dark 
matter candidates. The vectorlike coupling to the Z-boson results in an unsuppressed 
annihilation cross section into fermion pairs for low momenta, and in a large elastic scat- 
tering cross-section with nucleons. The latter strongly constrains the relic density of Dirac 
neutrinos to be several orders of magnitude smaller than the observed dark matter density. 

Some of these problems can be avoided when the Dirac pair is split using a Major ana 
mass term. In that case the axial coupling of A'^i to the Z boson helps avoiding most 
constraints from direct detection, provided that the mass splitting is larger than the typical 
0(100 keV) momentum transfer in dark matter nucleon scatterings. See e.g. ^30j for a 
discussion in the context of inelastic dark matter. 



pp e+re+r + ^ . 



(32) 
(33) 




(34) 
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5.1 Majorana Neutrino Annihilation to Fermions 



The annihilation rate of a Majorana fermion into standard model fermion pairs through 
^-boson exchange is well known from the case of neutralino annihilation in the MSSM. 
Here we briefly review the most important results. The coupling of A'^i to the Z-boson in 
four component language is proportional to 7^75. We consider the annihilation into a pair 
of SM fermions with mass nif and with a coupling to the Z boson of the form 

(T37m(1-75)-2Q/47m) , (35) 



where Qf is the fermion electric charge and Ts is the isospin eigenvalue. We work in the 

center of mass frame, where the A''i three-momentum can be replaced by |p| = y^s/i — . 
The spin-averaged squared matrix element, in unitary gauge, and integrated over solid 
angles, is then given by: 

dni^. . .,0 2cU^ 1 r3 



/ 



i:\M\^- ^ (,-M|;,.,|r| ^ [ i^f-K^ - Mlf. (36) 



where we used that Ts = ±1/2. This result displays the important helicity suppression for 
this process: In the |p| limit the amplitude is suppressed by the small fermions mass 
nif coming from the first term, while the second term vanishes as |pp. Further note that 
the momentum independent part is not resonantly enhanced at s = since the pole of 
the Z propagator is cancelled by the numerator. 

Prom this the annihilation cross section is obtained using 

A 2\ 1/2 

^rn}\ fdni,, 2 



M = — (37) 

Freeze-out typically occurs at temperatures T ~ Mi/20, such that the majority of Ni 
are non-relativistic. In this limit the center of mass energy can be approximated by s = 
4M^ -{-v'^Mi, where v is the relative velocity between the annihilating A'^i. Away from the 
Z pole region we can further set Tz — 0. Decomposing the annihilation cross section as 
(va) = a + bv^ , and dropping terms of order m^f^ and u^, we obtain the coefficients 

1 1 



1287r ci M| 

4g' 



Y^Nfm), (38) 



1 dg^ M2 / ^ 2 160 4 \ 



13 




50 100 150 200 

Mi[GeV] 



Figure 4: A^i relic density OjVi/i^ in the M1-M2 plane, for mh = 150 GeV and = 
Ml + 200 GeV. For comparison, the observed dark matter density is ^l^Mh'^ ~ 0.11. 



valid for Mi ^ ruf, and with Nf denoting the internal degrees of freedom of the fermion 
/. The result for b includes the sum over three families of quarks and leptons, excluding 
the top quark. The a term is suppressed by ml/M"^ relative to the b term, so that even at 
temperatures around the freeze-out the b term will dominate the annihilation. 
A good first estimate for the relic density can be obtained using 



1.04 X 10^ GeV"^ XF 



1 



where the freeze-out temperature xp = Mi/Tp is determined numerically from 



(40) 



xp = log 



c(c + 2 



45 Qd MiMpi(a + 66/xf) 




iQi^JXp 



(41) 



For our case we have gd = 2, g^, = 92 and we use Mpi = 1.22 x 10 GeV. The resulting relic 
density reaches at most 0(10%) of the observed dark matter abundance in the universe. 
For Ml ~ Mz/2 and Mi ~ m/i/2 the relic density is strongly suppressed by resonant 
annihilation into Z bosons or Higgs bosons respectively, while for Mi > Mw annihilation 
into gauge boson pairs becomes dominant and further reduces the relic density. 

In fig. [4] we show the A^i relic density in the Mi — M2 plane. To capture all possible 
annihilation channels and to include coannihilation effects, the model was implemented 
in micrOMEGAs 2.4 |32j . Clearly visible are the regions of resonant annihilation near 
Ml = 45 GeV and Mi = 75 GeV (the Higgs mass was set to m/j = 150 GeV). Furthermore 
the relic density is reduced close to the M2 = Mi line due to coannihilation, and for 
increasing Mi where annihilation into gauge boson and top quark pairs becomes possible. 
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5.2 Direct Detection of Majorana Neutrino Dark Matter 

The stringent direct detection constraints on Dirac neutrino dark matter are evaded pro- 
vided that M2 — Ml > 100 keV [30]. Furthermore if we assume that the A'^i abundance 
is purely thermal, then the direct detection rate has to be multiplied by the fraction of 
dark matter that is of A'^i type. In this case we find that the latest CDMS and Xenon 100 
results do not impose a constraint on the parameter space. 

In addition one can consider the case where all dark matter consists of A'^i. This can 
be achieved by reducing the coupling of A^i to the Z-boson, as suggested in [33]. Here we 
do not consider this possibility since this also affects the direct detection rates. Another 
way to increase the relic density of A'^i is to invoke a nonthermal production mechanism, 
e.g. using a heavy scalar that couples to the standard model only through a very small 
coupling to A'^i, and decays to A'^i pairs after freeze-out. It requires a somewhat accurate 
adjustment of masses and couplings to ensure that one obtains the correct relic density 
and an acceptable velocity distribution for A^i. 

Assuming that such a mechanism is at work, the elastic A'^i-nucleon cross section is 
again evaluated using micrOMEGAs. Both the spin independent (SI) and spin dependent 
(DS) scattering cross sections with nucleons turn out to be constrained by experiments. 

The currently most stringent constraint on SI nucleon scattering comes from the CDMS 
|34j and Xenon 100 |35j experiments, and have reached a sensitivity of 4 x 10^^ pb. The 
SI cross section is dominated by Higgs boson exchange and thus depends on m/i and on 
the light quark form factors of the nucleon. In the MSSM the strange quark form factor 
fs is a major source of uncertainty, since the couplings to down type quarks are enhanced 
by tan /3, see e.g. [36l [371 EB] for recent discussions of this issue. 

For the case considered here, the strange quark form factor is not as dominant, but 
still affects the exclusion limit, as will be shown later. The parameter dependence of the 
Higgs mediated scattering cross section is easy to understand. The A^i couplings to the 
Higgs boson are proportional to Mi and also enhanced for large M2, while the exchange of 
the Higgs boson in the t-channel gives rise to a mj^^ behavior. For small ruh this leads to 
exclusion of Mi values above 150 — 200 GeV. In fig. [5] the excluded region in the Mi — M2 
plane is shown for the case of a 150 GeV Higgs boson. The dependence on the strange 
quark form factor is illustrated by showing the constraint for three different values of fs- 

Constraints on spin-dependent neutrino-nucleon scattering are much weaker when com- 
pared naively to the SI constraints, the best limit coming from the Xenon 10 experiment 
[39j . It constraints the SD neutron cross section to be smaller than 5 x 10~^ pb. This 
limit turns out to be very restrictive for Majorana neutrino dark matter, due to its large 
coupling to the Z-boson. A purely Majorana neutrino is excluded in the mass range of 
10 GeV< M < 2 TeV. 

In the mixed Dirac-Majorana case considered here, the coupling to the Z-boson is 
reduced by Cg, thus the cross section is reduced by > j. In the case of maximal mixing, 
the constraint becomes somewhat weaker, but still excludes Majorana neutrino dark matter 
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Figure 5: Direct detection constraints on A'^i dark matter in the M1-M2 plane, for nih = 
150 GeV and = Mi + 200 GeV, under the assumption that = 0,bm- The red 
shaded areas are excluded by constraints on the spin dependent nucleon cross section, for 
the indicated values of fs- The whole parameter space is also excluded by constraints on 
the spin dependent neutrino-neutron cross section, as indicated by the blue stripes. 



for masses roughly between 15 GeV and 500 GeV, and thus the whole parameter range in 

In conclusion, forth generation neutrino dark matter is heavily constrained by dark 
matter direct detection experiments, even in the mixed Dirac-Majorana scenario, under 
the assumption that the neutrino relic density matches the observed dark matter abun- 
dance. These constraints are relaxed if either the local dark matter density and velocity 
distribution is different from the one assumed to obtain the limits, or if the Majorana 
neutrinos are only one of several dark matter species. In particular case where A^i is a 
thermal relic is not constrained by spin dependent neutron scattering, since the density is 
at least a factor of 5 smaller. This regime will however be probed by the next generation of 
experiments [40j. Finally A'^i annihilation in the sun and in the galactic halo is suppressed 
by the small thermal relic abundance, and therefore does not impose further constraints 
on the scenarios considered here [4:1\ [32] . 

6 Conclusions 

We have analyzed the effects of stable fourth generation Majorana neutrinos on Higgs 
phenomenology and on cosmology. 

The rate ahww = cr{pp — )• /i) x Br(/i — )• ly+VF") is an important test of the standard 
model of physics, and will eventually be measured at the LHC, provided that the Higgs 
boson is there. A deviation of the SM prediction is a clear indication that new physics 
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is changing either the Higgs production cross section or its branching fractions, or both. 
We have shown that in the presence of stable Majorana neutrinos, this rate (Jhww can be 
reduced significantly, despite the strong enhancement of the Higgs production cross section 
due to the fourth generation quarks. 

The semi-invisible decays of the Higgs boson into N1N2 and N2N2 pairs provide unique 
signals at the LHC, and might play an important role in understanding deviations in 
ahww- When the mass splittings between the neutrino states are small, these signals are 
very challenging, since the produced leptons are rather soft. 

A stable A^i is not excluded by current dark matter direct detection experiments, pro- 
vided that M2 — Ml is large enough to prevent large inelastic scattering rates. In this case 
A'^i will have a nonzero thermal relic density that can reach up to 20% of the observed 
dar matter density, but is below 10% for most of the parameter space. The possibility 
that all dark matter is of A^i type, e.g. invoking a non-thermal production mechanism, is 
excluded by constraints on the spin dependent dark matter-neutron cross section. Future 
direct detection experiments will strengthen these bounds and constrain the fraction of 
dark matter that can be of A^i type to be less than 10%. 

A stable fourth generation has an abundance of interesting phenomenological conse- 
quences. Besides the effects on Higgs phenomenology presented here, there can be signals 
of stable charged particles [15] , quark-antiquark bound states [ISj HH US] and lepton flavor 
violation [36] at the LHC. In the context of the present discussion, it would be good to 
investigate further the prospects of the LHC to detect the semi-invisible decay modes of 
the Higgs boson. 
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A Neutrino Couplings in Four Component Formalism 

A 2-component chiral field is stacked up with its Hermitian conjugate to form a 
4-component Majorana field = N (B ■ 
Thus we have 

+ h.c. Nl^lN^^ , (42) 

A^l + h.c. N^^IN^ , (43) 

iN2Ni + h.c. Ni^-f^Nf^ = NfU-f^Ni^ . (44) 
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As usual ip = Tp^'y^. Using the above relations it is easy to translate the Higgs boson 
couplings (10) into four component language 

£ = I'^h {2ceseNfNf + 2ceseNi' N^' + 2{sl - cDNfUj.N,^') , (45) 

The four component couplings to gauge bosons are obtained in a similar manner. For this 
we note that, in the Weyl basis, the 4x4 gamma matrices are defined as 



such that the lefthanded two component coupling Nja^Nj can be rewritten as 

iV/a^iVj N^^^j'^PLNf'^ . (47) 



For i = j only the axial part of the coupling survives. Therefore the couplings to the 
Z boson translate to 



9 1 r 



-ciiVf 7^f< - slN^'y^^N^^ + seceiNf'y^N^'\ . (48) 



Further useful relations between the masses M\^i and the quantities ce, sg are collected 
here: 

sin 6* m4Z) Mi [M[ 

tant' 



cos 6* M2 V -^2 



Ml / Ml m4D / M2 

sm f = — = = \ / — — — , cos ( 



Ml^m\^ VM1 + M2' JMf + mL V Mi + M^ ' 



. \/MiM2 2. .2^ M2-M1 

sm d cos = — -— , cos {) — sm " — 



Ml + M2 Ml + M2 

Finally the triangle function A that appears in the two body phase space integrals is defined 
as 
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